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ABSTRACT 

Order of magnitude estimates of radiogenic heat and antineutrino production 
are given, using the San Marco cathedral as an example. Prospects of determin- 
ing the radiogenic contribution to terrestrial heat by detection of antineutrinos 
from natural radioactivity (geoneutrinos) are discussed. A three kton scintillator 
detector in three years can clearly discriminate among different models of terres- 
trial heat production. In addition, the study of geoneutrinos offers a possibility 
of improving the determination of neutrino mass and mixing, by exploiting the 
knowledge of Th/U abundance in the Earth. 

1. A few facts about San Marco 

After the visit of the San Marco cathedral, accompanied by a most interesting 
historical and artistic overview, nothing should be added but our gratitude to Milla 
Baldo Ceolin for organizing this most interesting and timely conference and for spicing 
it with exceptional events. Physicists however are tenacious, and we cannot resist 
adding a few additional information, some of these a touristic guide will never tell 
you. 

First of all, San Marco contains radioactive materials. The cathedral mass being 
in the range of 100 kton , we expect it contains about 100 kg of Uranium, as the typical 
Uranium abundance in rocks is in the range of one part per million, however with 
large variations. We also expect some 400 kg of Thorium, since almost everywhere 



in the solar system (meteorites, Moon, Venus and also Earth) the typical abundance 
ratio is Th/U ~ 4. In addition, there are about 100 Kg of 40 K, corresponding to the 
typical ratio K/U ~ 10 4 which one finds in Earth rocks and to the natural abundance 
m K/K= 1.2- 10~ 4 . 

San Marco is also a heat source. In fact, each decay chain releases energy over 
long time scales {e.g. for 238 JJ A = 52 MeV and T\/2 = 4.5 Gyr ), the total heat flow 
being: 

H = 9.5M(U) + 2.7M(Th) + 3.QM( m K) (1) 

with masses in 100 kg and H in mW. This gives 24 mW for San Marco, really a very 
weak heat source, although it does not matter in these sunny days. 

More interesting to people attending this conference, San Marco is an anti-neutrino 
source. Each decay chain releases antineutrinos together with heat, with a well fixed 
ratio, (e.g. 238 JJ _> 206 pjj _|_ S A He + 6z/ + 52 MeV). The antineutrino luminosity is: 

L = 7AM(U) + l.QM(Th) + 27M{ m K) (2) 

where again M is in 100 Kg and L in 10 9 u/s. This gives about 4 ■ 10 10 z//s for San 
Marco. Antineutrinos from the progenies of Uranium (E max = 3.3MeV) and Thorium 
(E max = 2.2MeV) can be detected by means of inverse beta decay reaction: 

v + p^n + e + - 1.804 MeV (3) 

At least in principle, the two components can be discriminated, due to the dif- 
ferent end points. Unfortunately, antineutrinos from (3 decay of m K are below the 
threshold for (J3J), whereas neutrinos from electron capture are obscured by the Sun. 
The cross section of <jBJ) for 2.5 MeV antineutrinos ( a ~ 5 ■ 10~ 44 cm 2 ) corresponds 
to an interaction length in water A = 7 ■ 10 18 m. There is plenty of water near the 
cathedral and San Marco square (S ~ 10 4 m 2 ) is often covered with Acqua Alta 
(high water), a 10 cm height corresponding roughly to 1 kton, the size of KamLAND. 
Should Acqua Alta reach the clock of Torre dell'Orologio, the size of Superkamiokande 
is reached. With a height of 100m a megaton detector is obtained, the cathedral be- 
ing now deeply submerged. A pointlike source emitting 10 10 u/s with E=2.5 MeV at 
the center of a megaton water sphere will produce one event every four years. Also 
in view of the environmental impact of such a project, better we look at some other 
direction. 

2. The sources of terrestrial heat 

Earth re-emits in space the radiation coming from the sun (K & = 1.4 kW/m 2 ) 
adding to it a a tiny flux of heat produced from its interior ($ ~ 80mW/m 2 ). By inte- 
grating this latter over the Earth surface one gets a flow He = 407W, the equivalent 
of some 10 4 nuclear power plants. 



The origins of terrestrial heat are not understood in quantitative terms. In 1980 
J. Verhoogen concluded a review of terrestrial heat sources by saying^: "...what 
emerges from this morass of fragmentary and uncertain data is that radioactivity itself 
could possibly account for at least 60 per cent if not 100 per cent of the Earth's heat 
output. If one adds the greater rate of radiogenic heat production in the past, possible 
release of gravitational energy (original heat, separation of the core...) tidal friction ... 
and possible meteoritic impact ... the total supply of energy may seem embarassingly 
large. . . ". One can appreciate the complexity of the problem by comparing Sun and 
Earth energy inventories. In fact, a constant heat flow H can be sustained by an 
energy source U for an age t provided that U > Ht. For sustaining the sun over an 
age t so i ar = 4.5 • 10 9 yr, gravitational and chemical energies are short by a factor 10 2 
and 10 6 respectively and only nuclear energy can succeed, as beautifully demonstrated 
by Gallium experiments in the nineties. On the other hand, the terrestrial heat can 
be sustained over geological times by any energy source, be it nuclear, gravitational 
or chemical. 

Observational data on the amounts of Uranium, Thorium and Potassium in Earth 
interior are rather limited, since only the crust and the upper part of the mantle 
are accessible to geochemical analysis. As U, Th and K are lithofile elements, they 
accumulate in the c ontin ental crust (CC). Estimates for the Uranium mass in the 
crust are in the rang 

M C (U) = (0.2 -0.4)10 17 % . (4) 

Concentrations in the mantle are much smaller, however the total amounts are 
comparable due to the much larger mass of the mantle. Estimates for the whole 
mantle are in the range^ 

M m (U) = (0.4-0.8)10 17 % . (5) 

One has to remark, however, that these estimates are much more uncertain than 
for the crust as they are obtained by: i)collecting data for upper mantle (h u = 600km), 
ii) extrapolating them to the completely unexplored lower mantle (hi = 3000km). 

Concerning the abundance ratios, one has generally Th/U ~ 4, consistent with 
the meteoritic value. A remarkable exception is the oceanic crust where Th/U ~ 2, 
however both U and Th abundances are an order of magnitude smaller with respect 
to CC, which is also much thicker. 

Concerning Potassium, generally one finds K/U ~ 10,000. Earth looks thus 
significantly impoverished in Potassium with respect to Carbonaceous Chond rites , 
and also to other meteorites. This has long been known as the Potassium problenP^. 
In fact, elements as heavy as Potassium should not have escaped from a planet as 
big as Earth. It has been suggested that at high pressure Potassium behaves as a 
metal and thus it could have been buried in the Earth core, where it could provide 



the energy source of the terrestrial magnetic field. However, Potassium depletion is 
also observed in Moon and Venus rocks. The most reasonable assumption is that it 
volatized in the formation of planetesimals from which Earth has accreted. 

In conclusion, the determination of the Uranium, Thorium and Potassium in the 
Earth is an important scientific problem, as it can fix the radiogenic contribution to 
terrestrial heat production. 

3. Antineutrinos from below 

u If there are more things in heaven and Earth than are dreamt of in our natural 
philosphy, it is partly because electromagnetic detection alone is inadequate" . With 
these words in 1984, Krauss, Glashow and Schramm® proposed a program of antineu- 
trino astronomy and geophysics, which could open vast new windows for exploration 
above us and below. Now that we understand the fate of neutrinos it is time to tackle 
the program (including Earth energetics, a detailed study of the solar core, neutrinos 
from past super novae...) Determination of the radiogenic contribution to terrestrial 
heat production is the first step. 

One can build several models for the radiogenic heat production. Since for any 
element there is a well fixed ratio heat/(anti)neutrinos each model also provides a 
prediction for the antineutrino luminosities, the basic equations being (JTJ) and (J2J), 
where the same numerical coefficients can be used when masses are in units of 10 17 
kg, powers are in TW and luminosities in 10 24 /s. At this level, everything is fixed 
in terms of three inputs. The range of plausible models is covered in figQ which 
deserves the following comments: 

i) In a simple chondritic model one assumes that Earth is obtained by assembling 
together the same material as we find in these meteorites, without loss of heavy 
enough elements. The amounts of U, Th and K are determined from their ratio to 
Si in meteorites, by rescaling to the known abundance of this latter element in the 
Earth. This simple model easily accounts for 3/4 of He , mainly supplied from 40 A, 
however it implies K/U = 7 ■ 10 4 , a factor seven larger than the value observed in 
Earth rocks. 

ii) In the standard model of geochemists, the so called Bulk Silicate Earth (BSE) 
model, with K/U = 1 ■ 10 4 , the radiogenic production is about one half of the terres- 
trial heat flow, being supplied mainly from U and Th. 

iii) A fully radiogenic model, where the BSE abundances are rescaled imposing 
by imposing H rad = 40 TW, is not excluded by observational data. 

In all models, see fig 12 antineutrino production is dominated by 40 A decays. Th 
and U anti-neutrino luminosities are in the range (10 — 20) ■ 10 24 /s. By dividing over 
the earth surface, one gets fluxes of order 10 6 cm _2 s _1 , in the same range as that of 
solar boron neutrinos. 

Clearly, in order to estimate fluxes at a specific site one needs assumptions about 
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Figure 1: Contributions to terrestrial heat production (TW), as estimated in different models of 
Earth interior. CON= chondritic, BSE=bulk silicate Earth, RAD= fully radiogenic. The non 
radiogenic contribution is indicated as NR. 

the distribution of the radioactive materials, the total amounts being an insufficient 
information. In a paper submitted on December 2002 we provided estimates corre- 
sponding to different models for several sites''''. The numbers of predicted events for 
Kamland (normalized to an exposure of 0.14- 10 32 p-yr, a detection efficiency e = 78% 
and a survival probability P ee = 0.55) are 3.5, 4 and 6 for the chondritic, BSE and 
fully radiogenic model respectively. 

The paper ended with: "The determination of the radiogenic component of ter- 
restrial heat is an important and so far unanswered question. ... , the first fruit we 
can get from neutrinos, and Kamland will get the firstlings verusoon". 

A few days later the first results from Kamland appeared^, containing a first 
glimpse of Earth interior in addition to important information on neutrino oscilla- 
tion. Out of a total of 32 counts, in the prompt energy region below 2.6 MeV, 20 
events are associated with antineutrinos from reactors and 3 correspond to the esti- 
mated background. The remaining 9 counts are the first indication of geo-neutrinos! 
Clearly the uncertainty is large, since the expected statistical fluctuations are about 
y/32. Within this error, the result is consistent with any model for the radiogenic 
contribution to terrestrial heat, from to 100 TW. One has to remind that data have 
been collected in just six months and that significant accumulation can be achieved 
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Figure 2: Neutrino luminosity. Antineutrinos production rates from U, Th and K according to 
different models of Earth interior. In the last column neutrinos from ex. of K are shown. Units 
are 10 24 s _1 . 

in the future. 

4. Prospects for measuring terrestrial heat with geoneutrinos 

The first Kamland results really open a new window on Earth interior. A natural 
question is thus: what can be learnt on heat production from geoneutrino measure- 
ments? 

Let us consider Uranium geoneutrinos, since expected event numbers are higher 
and energies can be chosen in a range where background is negligible. The basic 
equations relating event number N(U), Uranium mass in the crust M C (U), in the 
mantle M m (U), and Uranium contribution to heat production rate H(U) are: 

H(U) = h[M c (U) + M m (U)} (6) 

N(U) = n c M c (U) + n m M m (U). (7) 

The first equation simply adds the contribution of crust and mantle, the propor- 
tionality constant being h = 9.5TW/10 17 kg. The different coefficients n c and n m 
reflect the different distances of the detector from sources in the crust and in the 
mantle. They also depend on the detector characteristics and on neutrino properties: 
as a normalization we shall consider an exposure of 10 32 p ■ yr a , a 100% detection 
efficiency and a neutrino survival probability P ee = 0.55. 

Assuming a uniform distribution in the mantle one has, independently of the site, 
n m = 18/(10 17 kg). On the other hand, n c depends on the location. Our preliminary 

a lkton of mineral oil contains 0.86 • 10 32 free protons. 




calculations, based on a crustal map of the whole Earth, give n c = 52/(10 17 kg) and 
n c = 63/(10 17 kg) for the Kamioka mine and the Gran Sasso respectively. 

There are two unknowns in eqs.(JBJ) and (J7jl . M C (U) and M m (U). From a single 
experiment one cannot extract both (in principle this could be achieved with two 
experiments, at markedly different geological sites, i.e., with very different values of 
n c ). 

One has to remind that M m (U) has been estimated by extrapolating to the whole 
mantle rather uncertain data obtained from the upper mantle. A reasonable program 
is thus: i)use available geochemical information to fix M C (U) and ii) use geo-neutrinos 
to get information on M m (U), and thus on H. In this way, from the above equations 
one gets: 

if (17) = (h/n m ) [N(U) - (n c - n m )M c (U)\ (8) 

A significant experiment should be capable of discriminating among models of 
heat production. In the BSE model one expects H(U) = 8.6 TW. Predictions for the 
Uranium contribution are between 6 and 16 TW, the lower (upper) value correspond- 
ing to the chondritic (fully radiogenic) model. A significant experiment should thus 
provide a measurement with a (la) error of about 2-3 TW. 

By considering the full range of estimated Uranium mass in the crust as a 3cr 
interval, one has M C {U) = (0.30 ± 0.07) 10 17 kg. From eq.(JBJ), which tells how to 
get H(U) once that N(U) has been measured with geoneutrinos, one derives the 
uncertainty on the extracted value of H(U) : 

AH(U) = (AH* + AHlf' 2 , (9) 

where: 

AH x = {h/n m )AN{U) , (10) 

and 

AH 2 = h {jlc ~ nm) AM C {U) . (11) 

Geoneutrino events are obtained from the counts number C after subtracting 
reactor R and background B events, so that: 

AN = C 1/2 = (N + R + B) 1/2 . (12) 

As demonstrated by KamLAND, background can be negligible in the energy region 
where most of U-events occur. The expectations for various experiments at different 
sites, collected in table [U suggest the following comments: 

i)KamLAND can reach an exposure corresponding to 10 32 p ■ yr (and full effi- 
ciency) in a rather short time. With this exposure, unless the nearby reactors power 
is reduced, the total uncertainty on H(U) would AH{U) ~ 8 TW, dominated by 
fluctuations of reactor events. A ten times large exposure (achievable in several 
years, with KamLAND, or in a much shorter time using a mineral oil detector with 



Superkamiokande size) becomes most interesting. One can estimate an accuracy of 
about 3 TW even with the present reactor power. This would really provide signif- 
icant information on the radiogenic heat production. Note however that some 270 
Uranium events will have to be extracted from a total of about 2000 counts and that 
in the present calculation accidental background has been neglected. 

ii) At Gran Sasso, where the reactor antineutrino flux is much smaller, an exposure 
corresponding to 10 32 p • yr and full efficiency should be reached in a reasonable time 
by Borexino. The uncertainty would be AH(U) ~ 5 TW. An experiment with ten 
times more data could reduce the uncertainty to about 2 TW. 

iii) As a limiting case, let us consider and ideal site where crust and mantle have 
the same weight in the determining neutrino events, so that N(U) is proportional 
to the total Uranium mass in the crust plus mantle. The ideal place is thus such 
that n c = n m , so that uncertainty on M C (U) do not affect the result, see eq.(JTTJ). 
Presumably this means a place in the middle of oceans, far away from the continental 
crust. This place is presumably also remote from nuclear plants, so that the error on 
N is essentially due to statistical fluctuations. Already with 10 32 p • yr one can reach 
an accuracy of about 2 TW. 

In conclusion, there are good prospects for reaching an accuracy AH(U) = (2 — 3) 
TW and thus fixing an important missing point on the sources of terrestrial heat. 

Table 1: The achievable accuracy AH on the determination of U contribution to heat flow. Calcula- 
tion are performed for 10 32 p- yr, 100% efficiency and a survival probability P ee = 0.55. Background 
counts are set to zero. In the last two lines, calculations are rescaled for an exposure of 10 33 p • yr, 
the achievable accuracy being now AH (33). Errors on heat flows arc in TW. 
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TV(best) 


26.4 


29.7 


16.2 


Geo-events from best BSE estimate 
(M C (U) = 0.3 and M m (U) = 0.6) 


C = R + B + TV(best) 


206.4 


64.7 


16.2 


Total counts 


AN = VC 


14.4 


8.04 
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Statistical fluctuations only 


AHx = (h/n m )AN 


7.6 


4.2 


2.1 


Error from geo-neutrinos 


AH 2 = h(n c - n m )/n m AM C 


1.3 


1.7 





Error from crust 


AH = y/AH'f + AH'i 


7.7 


4.6 


2.1 


Total error 


AHi(33) = (h/n m )AN 


2.4 


1.3 


0.7 


for lO^p ■ yr 


AH(33) = ^/AHf + AH'i 


2.7 


2.1 


0.7 


for 10 33 p ■ yr 



5. Can we learn on neutrinos from geoneutrinos? 



Uncertainties on the individual Uranium and Thorium fluxes $ are large. In prac- 
tice, one cannot use <&(U) and/or $(T/i) as additional input for extracting neutrino 
mass and mixing from the lowest energy region (E vis < E geo = 2.6 MeV) in a reactor 
experiment. 

However, one can still gain some information on 9 and Am 2 by observing that the 
ratio of events from U and Th is well constrained: 

r = N(Th)/N(U) = 0.25 ±0.05 . (13) 

This follows from the fact that the abundance ratio Th/U is well fixed in the solar 
system, being very similar on meteorites, Venus, Moon as well as on Earth , all this 
information pointing to a common origin of the solar system. 

The 20% uncertainty corresponds essentially to uncertainties on the Th/U ratio in 
the various components of the Earth which contribute appreciably to neutrino produc- 
tion (continental crust and mantle). This uncertainty accounts for the comparison of 
different estimates and also of different regions of the Earth interior. With respect to 
the average value Th/U = 3.8 value, the extrema are Th/U = 2 for the oceanic crust 
(which is however poor in both U and Th and also it is thin, so that it contributes 
little to neutrino production) and Th/U ~ 6 for the crust estimate by one author, 
most authors giving values in the range 3.8-4.2. An estimate Th/U = 3.8 ±0.7 should 
be more accurate than la. 

The constraint in eq.(|13jl was derived in^^ assuming a uniform Th/U distribution 
inside the Earth and assuming that the survival probability of geo-neutrinos reaching 
the detector is the same, independently of energy and distance. However, as remarked 
in^, the constraint holds within its error even if these simplifying assumptions are 
relaxed. 

Concerning the effect of regional Th/U variations, from^J it was found that r 
is changed by less than 2% when the detector is placed at Kamioka, or Gran Sasso, 
or Tibet (on the top of a very thick continental crust) or at the Hawaii (sitting on 
the middle of a thin, U- and Th- poor oceanic crust). Coming to the effect of local 
variations, by assuming that within 100 km from the detector the Uranium abundance 
is double, Th/U— 2 , one gets r = 0.22, whereas if its is halved, Th/U— 8, one finds 
r = 0.28. Neutrino oscillations clearly do not affect eq. (7) if the oscillation lengths 
for both U and Th neutrinos are both very short or very long in comparison with 
some typical Earth dimension. By numerical calculation, one finds that the effect of 
finite oscillation lengths does not change r by more than 2% for Am 2 > 1 ■ 10~ 5 eV 2 . 
In conclusion, all these effects are well within the estimated 20% uncertainty on r. 

The constraint (|T3*|) has been used in^^ in order to extract information from the 
full data set of KamLAND. A slight preference for the LMA-I solution with respect 
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Figure 3: The effect of the Th/U constraint on the determination of mass and mixing parameters, 
from ref. 10 

to LMA-II has been found, however it is not statistically significant. On the other 
hand, a significant reduction of the mixing parameters space has been found, in spite 
of the limited statistics. This constraint could become really useful when more data 
are available. 

6. Concluding remarks: from neutrons to neutrinos 

We all owe very much to Bruno Pontecorvo, the man who addressed the de- 
tectability of neutrinos, discussed available sources (Sun, reactors and accelerators), 
conceived the Cl-Ar method and invented the beautiful phenomenon of neutrino os- 
cillations. There is an additional lesson we can learn from him. In 1941, a few years 
after the celebrated Rome studies on slow neutrons, Bruno was a research physicist 
at Well Surveys Inc., and published a paper entitled "Neutron Well Logging - A New 
Geological Method Based on Nuclear Physics"^^. He had invented the neutron well 
log, an instrument still used for the prospection of water and hydrocarbons, see the 
page of the Society of Professional Well Log Analysts at http://www.spwla.org/. It 
consists of a neutron source and neutron or gamma detector (well shielded from the 
rays coming directly from the source) to be placed in the well. As hydrogen atoms 
are by far the most effective in the slowing down of neutrons, the detected radiation 



is primarily determined by the hydrogen concentration, i.e. water and hydrocarbons. 
In this way, the discoveries of the Rome group were applied to the study of quite 
different problems. Possibly now we have similar opportunities with neutrinos, due 
to the important achievements of the last few years: 

i) The fate of neutrinos is now essentially understood. We don't need anymore 
to rely on standard solar model calculations or on the comparison among different 
experiments. SNO has directly observed the transmutation of solar neutrinos and 
KamLAND has neatly confirmed the result with man made antineutrinos from nuclear 
reactors. 

ii) KamLAND has demonstrated that it can reach purity levels such that detection 
of geo-neutrinos is feasible, a most impressive experimental result. 

All this opens the road for using low energy neutrinos as real probes of nature, 
be it Earth, Sun, future and past supernovae as well as - may be - the big bang. 
Measurement of the radiogenic component of terrestrial heat will be the first step. 
It can be achieved with few kiloton (mineral oil) detectors in a few years, best if far 
from nuclear reactors. 

This goal can be reached without submerging San Marco inside megaton detectors 
and it will provide a definite answer to an important and long standing scientific 
problem. 
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